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ABSTRACT 

There is evidence that cosmic rays were present in galaxies at moderately high 
redshift. This suggests that magnetic fields were also present. If cosmic rays and 
magnetic fields must always be close to equipartition, as they are to an order of 
magnitude within the local universe, this would provide a powerful constraint on 
theories of the origin and evolution of magnetic fields in galaxies. We evaluate the 
role of magnetic fieldstrength in cosmic ray acceleration and confinement. We find 
that the properties of small scale hydromagnetic turbulence are fundamentally 
changed in the presence of cosmic rays. As a result, magnetic fields several 
orders of magnitude weaker than present galactic fields can accelerate and retain 
a population of relativistic cosmic rays, provided that the fields are coherent over 
lengthscales greater than a cosmic ray gyroradius. 

Subject headings: cosmic rays — galactic evolution — magnetic fields — MHD tur- 
bulence 



1. Introduction 

Superthermal particles, or cosmic rays, are a major constituent of the interstellar medium 
in galaxies. They are dynamically coupled to the thermal gas, and collectively provide pres- 
sure support, can drive outflows, and are thought to modify the structure of shocks. Cosmic 
rays are the dominant source of ionization, and a major source of heating, in regions of high 
visual extinction. Spallation reactions between cosmic rays and ambient interstellar nuclei 
are the primary mechanism for synthesizing light elements. 

The energy spectrum, directional anisotropy, and chemical composition of cosmic rays 
are measured directly through in situ observations on the Earth and in the heliosphere. The 
global properties of cosmic rays in the Galaxy are probed through their radio frequency 
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synchrotron emission and 7-ray line and continuum emission. A recent summary of the 
observations is given by Schlickeiser (2002). 

Although synchrotron and 7-ray radiation from normal galaxies have up to now been 
directly detected only within the local universe, there is evidence that cosmic rays existed in 
galaxies even at early times. The light elements present in metal poor halo stars are thought 
to originate by spallation reactions between energetic particles and the material from which 
these stars formed [Duncan, Lambert, & Lemke (1992), Gilmore et al. (1992), Duncan et 
al. (1998), Boesgaard et al. (1999) and references therein]. Further evidence is provided by 
the spectrum of the diffuse 7-ray background, which is best fit by a superposition of active 
and normal galaxy 7-ray spectra at moderate redshift [Vasiliki & Fields (2002)]. 

The energy densities in the Galactic magnetic field and cosmic rays are now roughly 
equal. If this relationship is universal and necessary, then the presence of cosmic rays in young 
galaxies implies the existence of equipartition magnetic fields. This would be an important 
constraint on theories of the history of magnetic fields in galaxies, a subject on which evidence 
is sparse. At present, the only direct observations of magnetic fields associated with young 
galaxies are the Faraday rotation measurements in damped Lya systems at redshift z < 1 
reported by Oren & Wolfe (1995). These observations suggest magnetic fields of ~ 1/xG 
strength which are spatially coherent over several kpc, similar to the fields in the Milky Way 
and other spiral galaxies. Inferences based on cosmic rays could probe galactic magnetic 
field evolution at earlier epochs, and, we will see, could be directly relevant to conditions in 
star forming regions. This is significant because of the important role magnetic fields can 
play in star formation and their possible effect on the Initial Mass Function. 

The leading theories of cosmic ray acceleration and confinement depend on the strength 
of the magnetic field. In this paper we use this dependence to estimate the minimum mag- 
netic field at which a population of relativistic cosmic rays with energy density comparable 
to the present Galactic population can be accelerated and confined. We find that a field 
several orders of magnitude weaker than the present Galactic field is sufficent for both accel- 
eration and confinement, provided that the coherence length of the field exceeds the cosmic 
ray gyroradius. 

The plan of this paper is as follows. In §2 we state our assumptions about those 
conditions in early galaxies which affect the acceleration and propagation of cosmic rays. 
In §3 we briefly review the theory of energetic particle transport and its application to 
cosmic rays. In §4 we derive constraints which arise by assuming that the cosmic rays diffuse 
at the minimum possible rate, which at a given fieldstrength gives the maximum acceleration 
and confinement. In §5 we discuss hydromagnetic turbulence in a weak magnetic field and 
find that it is drastically affected by even a small population of cosmic rays. Fluctuations 
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grow rapidly in the presence of cosmic ray anisotropy, and, in contrast to fluctuations driven 
by cosmic ray anisotropy in the present Galactic environment, are not Alfvenic in character. 
At small amplitudes, the damping rate is much lower than the excitation rate. Therefore, 
we predict that the waves grow to nonlinear amplitudes. This may have implications for the 
nature of turbulence in other environments in which the magnetic fields is relatively weak 
and a population of energetic particles is present. Section 6 is a summary and conclusion. 
Most of the technical material is in Appendices A - C. 

We use Gaussian cgs units throughout, except in expressing the energies of cosmic rays, 
where we follow standard practice and use eV, and also in certain formulae, in which we use 
"natural" units such as parsecs or years. 

2. Conditions in Early Galaxies 

We assume that the density and temperature of interstellar gas in young galaxies are in 
the range found in contemporary galaxies. At lower metallicity the cooling rate is reduced, 
and the period in which supernova remnants are nonradiative is extended. This moderately 
enhances the efficiency of cosmic ray acceleration. We assume that Type II supernovae are 
clustered, as they are now, and may occur at a higher rate than at present. We consider the 
overall galactic size scale to be similar to present galaxy sizes, or perhaps somewhat smaller. 
The values we choose are appropriate to spheroidal systems with radii comparable to the 
present thickness of galactic disks. 

The most important assumptions concern the magnetic field. A recent review of the 
origin of cosmological magnetic fields is given by Widrow (2002), and discussions of mag- 
netic fields in galaxies are given in Beck et al. (1996) and Zweibel & Heiles (1997). Briefly, 
theories of the origin of galactic magnetic fields can be classified as top down or bottom 
up. In top down theories, the fields arise through large scale coherent processes and have 
large scale structure, although they are generally very weak. Examples of top down the- 
ories include magnetogenesis during inflation (Turner & Widrow 1988), by the operation 
of the Biermann battery in cosmological shock fronts or other vortex structures [Pudritz & 
Silk (1989), Kulsrud et al. (1997), Davies & Widrow (2000)], in cosmological ionization 
fronts [Subramanian, Narasimha, & Chitre (1994), Gnedin, Ferrara, & Zweibel (2000)], 
and through the manyfold expansion of plasma lobes associated with active galactic nuclei 
[Furlanetto & Loeb (2001), Kronberg et al. (2001)]. Although only the first of these theo- 
ries produces fields which are truly coherent over cosmological scales, all three produce fields 
which have some coherence over galactic scales. 
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In bottom up theories, there are many independent, small scale sources of field. The 
magnetized plasma from these sources expands and diffuses to fill the interstellar medium, 
eventually merging to create pervasive fields with largely random structure. Examples of the 
sources include plerion supernova remnants and jets or winds from accretion disks surround- 
ing stellar mass compact objects [Zeldovich, Ruzmaikin, & Sokoloff (1983), Rees (1987)]. 
The magnetic field is inhomogeneously distributed until it has had time to diffuse. For ex- 
ample, if magnetic fields were seeded by massive stars, there could have been fields in OB 
associations before there was a large scale galactic field. 

Both top down and bottom up theories predict magnetic fields which are much weaker 
than present galactic fields, and must be amplified by several orders of magnitude or more 
to resemble the fields of several /iG which are observed in galaxies. It is thought that 
amplification is due to the action of a hydromagnetic dynamo, which stretches the field by 
fluid motions while dissipating it at the resistive scale. Since resistive effects typically operate 
at an AU or less, while the coherence lengths of galactic magnetic fields are observed to be 
several kpc or more, a longstanding issue is how the fieldlines can be lengthened, within a 
fixed volume, without the extensive folding and tangling that characterizes a predominantly 
small scale field [see e.g. Kulsrud & Anderson (1992) and Schekochihin et al. (2002)] 

For our purposes, the chief distinctions between top down and bottom up theories are 
that the former predict a large scale, spatially homogeneous field while the latter predict a 
small scale, inhomogeneous field. Partly in recognition of bottom up theories, we will discuss 
cosmic ray acceleration by supernova shocks and propagation in the general galactic field 
more or less independently of one another. 



3. Transport Theory with Application to Galactic Cosmic Rays 

There is a well developed theory for interactions between energetic particles and low 
frequency electromagnetic fluctuations Bi, E\ superimposed on a mean magnetic field -Bo- 
Reviews are given by Jokipii (1971), Wentzel (1974), Toptygin (1985), and Schlickeiser 
(2002). This theory has been reasonably successful in explaining the ~ 10 7 yr confinement 
time of cosmic rays, their near isotropy, and their energy density, the latter provided that 
approximately 10% of supernova energy is converted to cosmic rays. In the interests of 
keeping the paper self contained, we provide a brief treatment here. 

A particle with velocity v and cyclotron frequency u cr = ZeB /mcy interacts primarily 
with the Fourier components of B\ and E\ which satisfy the cyclotron resonance condition 

UJ — fc||1>|| ± UJ cr = 0, (1) 
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where u and k are the frequency and wavenumber of the Fourier mode and the subscript 
"||" denotes the projection along B . 

If the spectrum of waves is incoherent, if a particle interacts with a wave for one gyrope- 
riod, and if the effect of any single interaction is small, then the the wave-particle interaction 
produces diffusion in momentum space. Although the particles diffuse in both the magnitude 
and direction of their momentum p, the diffusion in angle is dominant as long as the phase 
velocity of the resonant waves is much less than the speed of light. The rate of angular 
scattering as a function of the variable p = p\\/p is 

Ka*) = \u cr {k R £(k^y l , (2) 

where S is the power spectrum of the magnetic field fluctuations, normalized to the large 
scale magnetic field, and kn is the value of k\\ which satisfies eqn. (1). 

The resonance condition eqn. (1) involves only the parallel component of k. However, 
if k± > k», the oscillations of the wave electromagnetic fields over the gyrorbit diminish the 
response of the particle, and the scattering frequency is reduced (Chandran 2000a). 

Locally, the fluctuations define a frame with velocity Uf, which is the mean fluctuation 
velocity measured in the rest frame of the observer and weighted by the intensity spectrum. If 
the fluctations propagate isotropically, or if the velocity u of the fluid is much larger than the 
relative velocities of the fluctuations and the fluid, then Uf can be approximated by u itself. 
If the scattering time v~ x and associated mean free path Ay = fiv jv are short compared to the 
global lengthscales and timescales of interest, the particles are well coupled to the fluctuation 
frame, with a small amount of diffusion caused by scattering. Under these conditions, the 
phase space distribution function f(x, p, t) of the cosmic rays is nearly isotropic with respect 
to p. The isotropic part of /, averaged over the ensemble of fluctuations, is governed by 
a transport equation which includes advection and compression by the fluctuations, spatial 
diffusion with diffusion tensor D, and a source S (x,p,t) 

(Skilling 1975). Particle drifts caused by large scale curvature or gradients in B can be 
included explicitly in eqn. (3), in which case they appear in the advection term (Forman, 
Jokipii, & Owens 1974). 

For our purposes, it suffices to consider a diagonal diffusion tensor with components D\\ 
and D±. The parallel diffusivity is related to the scattering frequency v by 

' 4 J-i <>{p) 
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According to eqns. (2) and (4), D\\ can be approximated by 



D \\ ~ \— {kR^ikR))' 1 ■ (5) 

6 tO r 



D ± ~ ; . „ J', ^-a - ( 6 ) 



The perpendicular diffusivity D±_ is given approximately in terms of D\\ by (Jokipii 
1987) 

D\\ 

l + (knJ(k R )y 

Equation (6) shows that cross field diffusion is weak if the fluctuation amplitude is small, 
but becomes nearly isotropic if the fluctuations are nonlinear {k R £(k R ) ~ 1). According to 
eqn. (2), if the fluctuation amplitude is nonlinear, the scattering frequency is approximately 
the cyclotron frequency. This is also the timescale on which a particle interacts with a single 
wave. 

Equation (3) is the foundation of the leading theories of cosmic ray acceleration and 
confinement. It gives the confinement time r c to a region of size L as 

T C ~^, (7) 

where D is the largest component of the diffusion tensor, generally D\\. Equation (7) is 
consistent with the observed confinement time of Galactic cosmic rays if k R £(k R ) ~ 1CT 6 — 
1CT 7 (Cesarsky 1980). Equation (7) is an upper limit in the sense that it does not include 
transport due to migration of the fieldlines themselves, which may be considerable if the 
lines of force have a large stochastic component (Berezinskii et al 1990). 

Equation (3) predicts that in a steady state in which diffusion of cosmic rays is balanced 
by sources, with little advection or compression, the relationship between the cosmic ray 
energy density U cr and the power density of the source V cr is given approximately by 

L 2 

U cr ~ r c V cr ~ P cr — . (8) 

In the Milky Way, U cr ~ 10~ 9 GeVcm -3 , and r c is known from the isotopic ratios of light 
elements to be ~ 10 7 yr. This implies that V cr ~ 5 x 10~~ 27 erg cm -3 s -1 , which is about 5% 
of the power density in supernovae. 

If the scattering is strong and the distribution function / is nearly isotropic, the cosmic 
rays can be treated as a fluid, and the dynamic and energetic coupling of the cosmic ray and 
thermal fluids, which is mediated by scattering, can be expressed in hydrodynamic form. 
The rate of momentum transfer F cr takes the form of a pressure gradient force 

F cr = -VP cr (9) 



-7- 



acting on the thermal fluid. The rate of energy transfer to the thermal fluid E cr can be 
written as 

E cr = -7cr(«/ ~ «) ' VP cr (10) 

(Volk, Drury, & McKenzie 1984), where 7 cr , the polytropic exponent of the cosmic rays, is 
usually taken to be 4/3. The rate of energy transfer is proportional to the frictional force 
between the cosmic rays and thermal fluid. The dynamical effects represented by eqns. (9) 
and (10) are responsible for hydrostatic support, acceleration of outflows, and modification 
of shock structure in the thermal interstellar gas. 

According to the leading theory of cosmic ray origin [original papers by Axford, Leer, 
& Skadron (1977), Bell (1978a), Bell (1978b), Blandford & Ostriker (1978); reviews by 
Blandford & Eichler (1987), Jones & Ellison (1991)], the particles are accelerated by the 
first order Fermi process operating in strong interstellar shocks driven by supernovae. Ac- 
celeration requires hydromagnetic turbulence both upstream and downstream of the shock. 
In a frame in which the shock is steady, the upstream turbulence approaches the shock at 
nearly the shock speed u s (the wave speed being negligible compared to the shock speed), 
and the downstream turbulence recedes from the shock at speed us/r, where r is the com- 
pression ratio. Particles which propagate upstream away from the shock, are reflected back 
across the shock, and are then scattered back upstream by the postshock turbulence gain an 
angle averaged net momentum of Am^usir — l)/3r, much as if they were trapped between 
converging walls. 

The steady state spectrum of particles accelerated by a parallel shock is a power law 
in momentum space, f(p) oc p~ q , where the spectral index q is a function only of r: q — 
3r/(r — 1). The power spectrum is realized only within a range p m in < P < Pmax- The 
lower limit p m in is the minimum momentum at which the conditions of shock acceleration 
theory are fulfilled: the particle speed v must exceed the shock speed us and Coulomb losses 
must be insignificant. The upper limit p max is the maximum momentum a particle can have 
reached since the onset of acceleration. In §4, we will use the expression for p ma x derived by 
Lagage & Cesarsky (1983) to derive a lower limit on the magnetic fieldstrength such that 
shocks can accelerate particles to relativistic energies. 

The original theory of shock acceleration predicts the shape of the energetic particle 
spectrum, but not its normalization. The latter depends on the mechanisms by which par- 
ticles are injected into the shock. This so-called injection problem is closely related to the 
issue of how the energetic particles modify the structure of the shock itself. We will parame- 
terize the injection rate by an efficiency factor representing the fraction of supernova energy 
which is converted to cosmic rays. Estimates based on theory and/or numerical experiment 
suggest that the efficiency is likely to be at least 10% [Ellison & Eichler (1984), Blandford 
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& Eichler (1987) Jones & Ellison (1991), Berezhko & Ellison (1999), Malkov, Diamond, & 
Volk (2000), Ellison (2002)] 



4. Constraints at the Minimum Diffusion Rate 

In this section, we assume that cosmic rays are scattered by nonlinear turbulence: 
k R £(k R ) ~ I. According to eqns. (4) and (6), the spatial diffusion coefficients D\\ and 
D± are approximately 

D\\~D±~ — 3.1 xlO 16 — -, (11) 

3uj cr B Z 

for particles with charge Z, atomic number A, speed (3c, and Lorentz factor 7. Equation 

(11) corresponds to scattering by one gyroradius r g per gyroorbit, which is generally thought 

to be the minimum rate at which cosmic rays can diffuse. We defer the question of whether 

such a high level of turbulence can be maintained to §5. Here, we merely find the minimum 

value of B which is consistent with the theories of cosmic ray acceleration and propagation 

at maximum efficiency. 

In this section and the following one, we focus on mildly relativistic cosmic rays with 
energies of about 1 GeV/nucleon. This is slightly more than the energy required for spallation 
reactions (a few hundred MeV/nucleon), and the energy required to generate n mesons from 
nuclear collisions (~280 MeV/nucleon). That is, we focus on the energy range required to 
create light elements and the 7-ray continuum resulting from tt° decay. We do not discuss 
the acceleration of electrons. They are underabundant in cosmic rays, probably because they 
are injected at lower efficiency, but the injection problem is beyond the scope of the present 
paper. 

The acceleration and propagation of the MeV cosmic rays responsible for ionization 
and heating the ISM is not, in our view, a plasma physics problem. The energy of these 
ions is comparable to the energy per nucleon of supernova ejecta, and their range against 
inelastic collisions in neutral hydrogen is 10 20 - 10 21 cm -2 . Therefore, these particles probably 
originate as supernova debris, and decay within their first traversal of the galaxy rather than 
being confined for times much longer than their transit time. Interstellar gas has probably 
been heated and ionized by low energy cosmic rays since the first generation of supernova 
explosions. 
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4.1. Acceleration 

Lagage & Cesarsky (1983) have derived the rate of particle acceleration by a strong 
parallel shock when the fluctuation amplitude is unity. Their result can be written as 

d_E = 3ZeB 2 = 10-^ZBulGeVs- 1 , (12) 
dt 20 c s s K J 

where B and us are expressed in G and cms -1 , respectively. 

Equation (12) can be integrated in time for any model of the shock velocity us(t). LC 
found that most of the energy gain occurs during the blast phase, when the remnant is 
essentially undecelerated and us is constant, and the energy conserving Sedov phase, when 
us oc t~ 3 / 5 . The duration t b of the blast phase is 

M 5/6 

tb = 6 ' 2 x 109 ^7^ S ' (13) 



where M e is the mass of the ejecta in M & , n is the ambient interstellar number density in 
cm -3 , and £51 is the energy released in the supernova explosion, in units of 10 51 erg. From 
eqn. (12), ^ 

E(t b ) = 9.1 x 10 9 ^ ZB GeV. (14) 

The energy gained between the beginning of the Sedov phase and its termination by radiative 
cooling at time t r can be written in terms of the shock radii R b and R r at times t b and t r 

AE S = 1.3 x 10 10 4C^" I 1 ~ (%) m ) GeV - ( 15 ) 

M l J\^ \ \Rr) ) 

If the metallicity of the ambient medium is low, R r will be somewhat larger than its contem- 
porary value, somewhat increasing the energy gain. In any case, Equations (14) and (15) 
show that in order to accelerate protons to relativistic energies, E ~ GeV, B must be at 
least ~ 1(T 10 G. 

Since Type II supernovae are observed to be clustered, we also consider a shock driven 
by expansion of a superbubble driven by a constant (radiative plus mechanical) luminosity 
L [Shull & Saken (1995)]. In this case, us oc t~ 2/5 and the bubble radius R s oc t 3/5 . The 
particle energy E sb increases steadily with time, and can be written in terms of the bubble 
radius Rs as 

E sb = 8.2 x 10 s {^^j ZB GeV, (16) 
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where L 38 and R pc are L in units of 10 38 erg s _1 and the bubble radius in pc, respectively. 
Equation (16) shows that particle acceleration by superbubbles is generally weaker than 
acceleration driven by single supernovae. A similar conclusion was reached by Bykov & 
Fleishman (1992) on the basis of a somewhat different acceleration model. 

If the upstream magnetic field lies nearly in the plane of the shock, instead of nearly 
perpendicular to it, particles drift along the shock face while scattering back and forth 
across it. The maximum energy which can be acquired by drifting a distance R is of 
order (Jokipii 1987) 

ZcB 

E d ~ — —u s R = 30ZBu s RpcGeV. (17) 

If us ~ 10 9 cms" 1 , and R ~ 1 pc, B must again be about 10~ 10 G in order to accelerate 
protons to 1 GeV. 

The requirement that the particle gyroradius must be less than the radius of curvature of 
the shock also imposes a lower limit on the magnetic fieldstrength. However, the gyroradius 
of a GeV proton in a 10~ 10 G magnetic field is ~0.01 pc, so the constraint arising from the 
acceleration rate is more severe. 



4.2. Confinement 

The constraint that the cosmic ray gyroradius must be less than the typical magnetic 
lengthscale L B in the galaxy is thought to set an upper bound on the energies of cosmic 
ray nuclei in contemporary galaxies, and has led to the inference that the highest energy 
cosmic rays are extragalactic in origin. This limit arises because if r g > L B , the cross field 
drift velocity becomes comparable to the particle velocity itself. The constraint derived from 
setting r g = L B is 

B>w -^Vp±, (18) 

where L B is expressed in kpc. 

The limit on B derived from eqn. (18) is well below the value derived from the accel- 
eration time [eqns. (14), (15), (16), or (17)], provided that L B is comparable to the scale 
of the galaxy itself. However, it is a nontrivial constraint on those dynamo theories which 
predict that much of the power in the magnetic field is on small scales, especially when the 
field is still very weak [Kulsrud & Anderson (1992), Schekochihin et al. (2002)]. If the 
power spectrum in the magnetic field peaked at the Ohmic scale, or even the viscous scale, 
a much stronger field would be required to confine cosmic rays. 
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If we require that cosmic rays are confined within superbubbles long enough to produce 
light elements by spallation reactions, the resulting condition on B is more stringent 2 . Using 
eqns. (11) and (7), the resulting constraint on B for particles with Lorentz factor 7 and speed 
v = [3c is 

B > 1.0 x 10- 7 !^G, (19) 

where tq is r c expressed in units of 10 6 yr. Equation (19) shows that a field of 1CT 11 G is 
sufficient to confine GeV cosmic rays for 10 6 yr in a 100 pc superbubble, provided that the 
turbulence is strong. This is comparable to, but slightly weaker than, the field required for 
acceleration of such cosmic rays. 

Equations (8) and (11) lead to an expression for B in terms of the energy density in 
cosmic rays U cr and the power density of the cosmic ray source V cr , which we write as a 
fraction e of the supernova power density V sn . Taking 7/3 2 = 1 as average, setting A = Z = 
1 because most cosmic rays are protons, normalizing U cr and P sn by their contemporary 
values: U cr = 10~ 9 w cr GeV cm~ 3 , V cr = 10 _22 p cr eGeV cm~ 3 s, and expressing the confinement 
lengthscale L in kpc gives 

B = 3.3x 10- 14 U 2 cr G. (20) 

eL kpcPcr 

Nonlinear models of cosmic ray acceleration in shocks suggest that e > 0.1. If all the 
other parameters take their fiducial values, a field of about 10~ 12 - 10~ 13 G would sustain 
a cosmic ray energy density comparable to that found in contemporary galaxies. An even 
lower fieldstrength would have sufficed if the supernova rate was much higher at these early 
times than it is now. 



5. Wave Propagation, Excitation, and Damping in Weak Fields 

In this section we consider the propagation, excitation, and damping of hydromag- 
netic fluctuations at very low magnetic fieldstrength, and attempt to estimate whether early 
galaxies could have sustained a nonlinear level of magnetic fluctuations. 

At the magnetic fieldstrengths considered here, the frequencies of Alfven waves at the 
resonant wavenumbers are extremely low. The Alfven speed vai defined with respect to 
the ionized mass density is vai = -B/(47rmjnj) 1 / 2 = 2.2 x 10 11 -B/n l 1 ^ 2 cm s _1 . Motivated by 



2 Unraveling the history of light elements from measurements of their relative abundances in the oldest 
stars is a complex matter which may require more than one type of spallation source; see Parizot & Drury 
(1999), Parizot & Drury (2000), Suzuki & Inouc (2002). We provide only a rough estimate here. 
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the resonance condition eqn. (1), we define a fundamental wavenumber ko in terms of the 
nonrelativistic proton gyrofrequency id cp 

fc = — = 3.5x 10~ 7 B cm- 1 . (21) 

c 

The Alfven wave frequency is 

B 2 

u ~ k v Ai = 7.3 x 10 4 — ^ s _1 . (22) 

rii 1 ' 1 

For example, if B ~ 10~ 10 G and rij ~ 1, the wave period is more than 10 8 yr. 

The implication of the very long Alfven wave period is that the waves are strongly 
affected by cosmic rays. Under contemporary interstellar conditions, cosmic rays influence 
hydromagnetic waves only through the resonant interaction, described by eqn. (1). This 
means that only the imaginary part of the wave frequency is affected by the cosmic rays. 
But when the magnetic field is weak, the electromagnetic response of the plasma is so feeble 
that even a low density of cosmic rays has a significant effect on both the real and imaginary 
parts of lj. A similar situation prevails if a very high flux of cosmic rays is present [Kulsrud 
& Zweibel (1975), Zweibel (1979)]. 

At low magnetic fieldstrength, Alfven waves are also affected by the thermal plasma. 
Foote & Kulsrud (1979) have shown that the the large thermal ion gyroradius causes the 
ions to deviate from the E x B drift which nearly exactly characterizes their motion in a 
cold plasma, leading to modifications to the Alven wave dispersion relation. Thermal plasma 
also causes efficient damping outside a small cone of propagation angles centered on Bq. It 
turns out that with the lower bounds on B imposed in the previous section, the thermal 
modifications to the real part of td are not as important as the cosmic ray modifications. 
We justify this statement below. Collisionless damping is potentially important, however, 
and we consider it in §5.2. Because the collisionless damping increases with the propagation 
angle 9 as 6> 2 , while the excitation rate due to cosmic ray streaming is relatively insensitive 
to 9 for small 9, we consider only excitation of waves with = 0. 

Much of the interstellar gas in early galaxies may have been only weakly ionized, leading 
to damping of the waves by ion-neutral friction. We consider weakly ionized gases in §5.3. 

In §2 we mentioned the possibility that young galaxies do not have large scale magnetic 
fields. In discussing wave propagation, we tacitly assume that there is a well ordered field on a 
scale at least several times larger than the wavelength of the resonant wave; A ~ 6 x 10 _12 i? _1 
pc. Thus, a field of 10~ 12 G would have to be coherent over tens of pc in order to support 
waves. This is much less than the current coherence length of the Galactic magnetic field, 
which is at least a few kpc. However, the possibility remains that even this degree of 
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coherence is lacking. In this case, it seems unlikely that the particles can be confined at all; 
see eqn. (18). 



5.1. Dispersion relation with cosmic rays 

We assume the cosmic rays are protons. As the presumptive dominant species, they 
control collective effects. The dispersion relation for hydromagnetic waves in a medium 
with cold protons and cosmic ray protons drifting relative to each other with speed vd, and 
with electrons cancelling the charge and current of each population of protons, is derived 
in Appendix A [eqn. (A10)]. The dispersion relation is modified by cosmic rays because of 
their drift and by the same mechanism that modifies the dispersion relation in a hot plasma: 
the electron motion in the wave is nearly the E x B drift, but the cosmic rays deviate from 
the E x B drift because of their large gyroradii. Therefore, the currents associated with 
perturbations of the two species do not come close to cancelling, and even at low density the 
cosmic rays have a strong effect on the waves. 

In the rest frame of the thermal plasma, the dispersion relation is 



00 



+ u cp — (ir(k) {oj - kv D ) - k v Ai = 0. (23) 



n 



The complex functions (i r {k) (the subscript denotes the sign of circular polarization) are 
defined in Appendix A and plotted for a sample cosmic ray distribution function in Figure 
1. The real and imaginary parts of ( represent the contributions of the nonresonant and 
resonant cosmic rays, respectively, and are of similar magnitude for cosmic rays near the 
mean momentum of the distribution. Because the (i r cannot in general be easily evaluated, 
we treat Q r as a parameter of order unity in what follows, except in computing the group 
velocity [see eqn. (39)]. If vd = vau oj is real, and the waves are marginally stable. 

It is convenient to scale the wavenumber by 

k = yk , (24) 

where ko is defined in eqn. (21), and to absorb the ratio of cosmic ray to thermal plasma 
density into the Q r by introducing 

Qr = Or — - (25) 
Hi 



With these definitions, the general solution of eqn. (23) is 



uo 1 

uj cp 2 



-e lr ±[ef r + + 4ei r y^y 



(26) 
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The properties of eqn. (26) depends on the parameter 

X = „ 1.5 x 10-10 ( l lj^L I ( — ) (27) 

V v A iV Al y yn\ /2 B " " 

which is the ratio of the third term under the radical to the second term. This parameter has 
a more physical interpretation: if we assume that most of the cosmic ray pressure is provided 
by mildly relativistic particles (7 ~ 1) and use eqns. (25) and (27), then the ratio of the 
pressure in the mean magnetic field to the cosmic ray pressure can be written approximately 

as 

8nP cr 21 c' 1 ; 

Equation (28) shows that if X is large, the magnetic field must be far below equipartition 
with the cosmic rays. 

Under contemporary interstellar conditions, X <C 1 unless the cosmic ray flux is much 
larger than average (Zweibel 1979), and eqn. (26) reduces to 

^ = ±y ^ + f( 1± ^). ( 29 ) 

UJ C p C 2 \ VAiJ 

Equation (29) agrees with the expression for the growth rate given in Kulsrud & Cesarsky 
(1971). 

The drift velocity vp, diffusion coefficient D, and cosmic ray gradient lengthscale L c are 
related by v D ~ D/L c . Setting D equal to is minimum value [eqn. (11)] gives the minimum 
possible value of vd/vai 

^~A^ 4 .7xlO- 17 ^, (30) 
^Ai Vvm L c kpc B 2 1 J 

where in the last step L c is expressed in kpc. At the fieldstrengths of interest here, B < 
1CT 10 G, the anisotropy produced by cosmic ray sources drives the instability far above thresh- 
old. Under these conditions, the solution of eqn. (26) is approximately 

u j-f v d\ 1/2 ,01 \ 

— = ±[e lr y — . (31) 

The real and imaginary parts of eqn. (31) are comparable in size, and the wave periods are 
much shorter than those implied by eqn. (22). Evaluating eqn. (31) numerically gives 



-c = 3.2 x KT 1 -^ (^) 1/2 (lOVvCir) 172 . (32) 
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Although the turbulence required for diffusive shock acceleration might always be present 
in the ambient medium, it is of interest to investigate the circumstances under which the 
streaming instability grows quickly enough for the cosmic rays to confine themselves. A 
rough estimate of n CT /rii follows from the efficiency assumed for shock acceleration and the 
ratio of the mean cosmic ray energy to the postshock thermal energy: if 10% of the shock 
energy is converted to cosmic rays and the energy ratio is 10 5 , then n cr /rii ~ 10~ 6 . If 
vd ~ us ~ c/30, eqn. (32) gives u ~ 1.9(y(i r ) 1 ^ 2 B. If B ~ 1CT 10 G, which is approximately 
the limit implied by eqns. (14), (15), and (17), the growth time of the waves is a few hundred 
years. Thus, it is not implausible that cosmic rays in the vicinity of a strong shock can trap 
themselves. 

In order to investigate confinement on the galactic scale, we express vc in terms of the 
confinement time T7, the confinement time r c in units of 10 7 yr, and Lj, pc , the confinement 
lengthscale in units of kpc. If q r ~ 1CT 9 , lut c ~ 10 12 5, which exceeds unity as long as 
B > 1CT 12 G. This is comparable to, although slightly larger than, the field required to 
maintain cosmic rays at their present energy density [eqn. (20)]. 

Finally, we return to the effect of thermal plasma on the wave frequency. According 
to the results of Foote & Kulsrud (1979), the ion gyroradius correction is represented by 
adding the term 

y 2 v 2 

±uu cp —-^ (33) 

to eqn. (23), where the ± signs again represent left and right circular polarization, respec- 
tively, and v 2 = 2k B T/rrii = 1.6 x 10 8 T. The relative effect of this thermal term on the 
dispersion relation is measured by the quantity 



y3 ^V^~3.4xlO- 26 T 2 (^M ( — ). (34) 



16e ir \ cJv D \n cr J \v D 

The right hand side of eqn. (34) is typically much less than unity, so we neglect thermal 
corrections to Re(uj). 



5.2. Landau damping 

Unless the angle 6 between the wave vector and B is exactly zero, the wave electric 
field has a component E\\ parallel to B . Thermal ions which satisfy the Landau resonance 
condition u = k\\v\\ are steadily accelerated by E\\ and absorb energy from the wave. This 
process is known as Landau damping, and is the primary linear damping mechanism in fully 
ionized regions. 
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As we argue below, the condition 9 = is unlikely to be fulfilled because of inhomo- 
geneity in the background density and magnetic field. Therefore, 9 ^ is the generic case, 
and we must consider damping. 

Foote & Kulsrud (1979) evaluated the rate of Landau damping in a hot thermal plasma, 
with vf/v\ { 3> 1. We show in Appendix B that their result is unchanged by the presence of 
cosmic rays. The damping rate 1^ is given to within a factor of order unity by 



Tl = ^fhvj 2 , (35) 

where we have approximated tan# by 9 [see eqn. (B4)] within the range of interest. 

We calculate the wave damping time in two ways. First, we assume that the waves are 
excited at 9 = 0, but 9 grows as the waves propagate through the inhomogeneous background. 
We define the Landau damping time by the condition 

rt L 

/ T L dt = 1, (36) 
Jo 

where t is the time along a ray path. 

Let L b denote the gradient lengthscale of the background, and v g the group velocity of 
the wave; v g = du/dk. The evolution equation for 9 is 

< 37 > 

with the solution 9 = v g t/L starting from zero initial condition. The solution of eqn. (36) is 



y/nk\\Vi V 2 



According to eqns. (31) and (25), 



( n cr u cp v D \ 1/2 

where xll 2 = R e [Ci r 1 ^ *d(kCir) I 'dk\. Using eqn. (39), we rewrite eqn. (38) as 



tL = * „ 3 Q x 1Q9 f^_±) " 3 s , (40) 

WnvDViXirVcpJ \ \u n cr "x,i r T I B vd 



where L h pc is L b expressed in pc. 



Damping and excitation balance when ujt^ ~ 1. Using eqns. (32) and (40) yields the 
condition 



Evaluating the right hand side of eqn. (41) with vd set equal to its minimum value given in 
eqn. (30) yields 



If u)ti\ m i n exceeds unity, Landau damping is always weaker than excitation by cosmic ray 
anisotropy. This is likely to be the case for B > 10 -11 — 10 _10 G. Presumably, the wave 
amplitude is limited to Bi/B ~ 1 by a nonlinear mechanism. 

As an alternative method of calculating the damping, we suppose that the initial con- 
dition 6 = can never be realized in an inhomogeneous system. If instead of solving eqn. 
(37) we simply set 9 = (k^L)^ 1 and define the damping time t' L as T^ 1 , with I\ taken from 
eqn. (35), we recover the condition on vd given by eqn. (41) with the factor of 1.1. reduced 
by about a factor of 2. The insensitivity of the damping rate to the initial conditions justi- 
fies our approximate treatment of the problem, and reinforces the conclusion that Landau 
damping by the thermal plasma is weak compared with excitation by cosmic rays. 



Large volumes of the interstellar gas in young galaxies may be only weakly ionized. In 
such regions, hydromagnetic waves are damped by ion-neutral friction. Under present con- 
ditions, frictional damping is so strong that cosmic ray stream through H I regions virtually 
without scattering (Kulsrud & Cesarsky 1971). 

The rate T in at which Alfven waves are damped by friction is given in Kulsrud & Pearce 
(1969). The modification of the waves by cosmic rays significantly reduces the damping rate. 
As shown in Appendix C, most of the energy in the waves is in the cosmic rays themselves, 
with a smaller fraction in electromagnetic fields, and the energy carried by the thermal gas 
a distant third. Since the ions are coupled to the wave by electromagnetic forces, while 
collisions with the neutrals directly tap only the small amount of energy carried by the ions, 
the overall effect of ion-neutral friction is doubly weakened. 

The damping rate is calculated by an energy method in Appendix C. The result can 
be written in terms of the ion-neutral collision frequency u in = p n (av) / \rrii + m n ) and the 




(41) 




(42) 



5.3. Ion-neutral damping 
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phase velocity = ui/k as 



We assumed in deriving eqn 
appears to be well justified. 

If the waves are Alfven waves, v<j, = vau and eqn. (43) reduces to v in /2. This is the 
usual expression for the damping rate of high frequency Alfven waves in a weakly ionized 
gas (Kulsrud & Pearce 1969). According to eqns. (27) and (31), 

^ ~ X 1 ' 2 » 1 (44) 

^Ai 

in the regime of interest. Not only is the frictional damping rate a factor of X 2 lower than 
it is for Alfven waves, it actually decreases as v D /c increases. Therefore, the final streaming 
rate cannot be determined by balancing ion-neutral friction against excitation by cosmic 
ray anisotropy, and the ionization state of the medium has little relevance to cosmic ray 
acceleration or confinement. 



(43) that uj > u in . In the parameter regime of interest, this 



6. Summary and Conclusions 

The light elements present in the oldest Galactic halo stars, and the spectrum of the dif- 
fuse 7-ray background, are indirect evidence that cosmic rays were present in young galaxies. 
Cosmic rays cannot be accelerated or confined without some level of magnetization. In this 
paper, we investigated just how strong a magnetic field is required, and what the presence 
of cosmic rays implies about the evolution of galactic magnetic fields. 

Cosmic rays can be confined in appreciable numbers only if the magnetic coherence 
length L B is larger than their gyroradii. This is the standard argument for the extragalactic 
origin of the cosmic rays with energies above ~ 10 9 GeV. Generalizing the argument to early 
times leads to a lower limit on the product BL B [eqn. (18)]. If L B were as large as a kpc, 
typical of the gradient lengthscale in the galaxy itself, then B must have exceeded 10~ 15 G 
if GeV protons were confined. If the magnetic energy density was concentrated at a much 
smaller scale such as the the Ohmic dissipation length, cosmic rays probably could not have 
been confined at all. 

Under contemporary conditions, cosmic rays are scattered by gyroresonant interactions 
with small amplitude (Bi/B ~ 1CT 3 — 1CT 4 ) hydromagnetic turbulence [eqn. (1)], and 
propagate diffusively. Although the turbulence can in principle be excited by the streaming 
anisotropy of the cosmic rays themselves, it is also strongly damped, and it appears that 
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supplemental sources of turbulence are required 3 . The excitation of waves by cosmic rays is, 
however, important in the vicinity of strong shock waves, where it is believed that the bulk 
of galactic cosmic rays are accelerated. 

If cosmic rays in early galaxies propagated diffusively, the limits on the fieldstrength 
are more stringent. We derived lower bounds on the fieldstrengths necessary for diffusive 
shock acceleration and for confinement by assuming that the scattering frequency is the 
gyrofrequency. This corresponds to a nonlinear level of turbulence [Bi/B Q ~ 1; eqn. (2)]. 
We used the upper bound on the acceleration rate derived by Lagage & Cesarsky (1983) to 
estimate the minimum magnetic fieldstrength required to accelerate protons to relativistic 
energies [eqns. (14), (15), (16), and (17)] in supernova and superbubble driven shocks, and 
found it to be about 1CT 10 G. This turns out to be similar to the minimum fieldstrength 
required to confine cosmic rays to superbubbles long enough for them to undergo nuclear 
collisions with ambient material and synthesize the light elements observed in the oldest 
stars [eqn. (19)]. At this fieldstrength, the growth rate of the streaming instability is fast 
enough that the cosmic rays may be able to trap themselves in the vicinity of the shock front 
[see discussion following eqn. (32)]. 

The magnetic fieldstrengths in regions of active star formation may not have been repre- 
sentative of the global galactic fieldstrength (see the discussion in §2). Therefore we derived 
a relationship between the average interstellar energy density in cosmic rays, the power in 
cosmic ray sources, and the global confinement time, assuming the maximum scattering rate. 
If the intensity of cosmic ray sources and the energy density in cosmic rays were comparable 
to their current values, the magnetic field could have been as low as 1CT 12 - 1CT 13 G. The 
fieldstrength scales inversely with the luminosity of the sources, so an enhanced supernova 
rate would have permitted an even smaller magnetic field. A field of this strength is consis- 
tent with growth of the waves on a timescale less than the nominal confinement time of 10 7 
yr, provided that the drift velocity v D exceeds ~100 km s -1 . 

Estimates of the fieldstrength based on cosmic ray diffusion theory scale with the tur- 
bulent amplitude as (Bx/Bq)^ 2 [see eqn. (2)]. There are two reasons why the turbulence 
could have been much stronger at early times than it is today. The sources extraneous to 
cosmic rays range from stellar radiative and mechanical luminosity to large scale dynamical 
effects associated with differential rotation and self gravity [Sellwood & Balbus (1999), Kim 
& Ostriker (2000)]. Under present conditions, the gyroradius of a GeV cosmic ray is about 



3 If the turbulence is highly anisotropic, so that the efficiency of scattering is reduced, then processes 
other than gyroresonant scattering, such as magnetic mirroring, may be necessary to explain confinement 
(Chandran 2000b) 
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1 AU [see eqn. (21)], much smaller than the scale at which turbulent energy is injected. But 
at fieldstrengths of 10~ 12 G, the gyroradii were at parsec scales, closer to the energy injection 
scale, at which the turbulent amplitude is significantly larger. 

The strength of these turbulent sources, however, can only be speculated upon. Tur- 
bulence is driven by the anisotropy of the cosmic rays themselves. The anisotropy is an 
inevitable consequence of the discrete and inhomogeneous distribution of cosmic ray sources, 
and the finite size of the galaxy. In §5 we considered the propagation, excitation, and damp- 
ing of low frequency electromagnetic waves in a medium consisting of cosmic rays, thermal 
plasma, and a weak magnetic field. We found that the circularly polarized, parallel prop- 
agating Alfven waves which efficiently scatter cosmic rays in contemporary galaxies, and 
which are excited by cosmic ray streaming anisotropy, are dramatically modified because 
the electrons E x B drift in the wave electromagnetic fields but the cosmic rays do not 
[Appendix A and eqn. (27)]. These modifications are only significant when the pressure in 
the mean field is much less than the cosmic ray pressure [see eqn. (27)]. We found that if the 
drift velocity vd is far above the instability threshold vau the real and imaginary parts of the 
wave frequency are comparable, and much larger than the Alfven frequency [see eqn. (31)]. 
In contrast to the current situation, neither Landau damping in fully ionized regions (§5.2 
and Appendix B) nor frictional damping in weakly ionized regions (§5.3 and Appendix C) 
appears capable of competing with the excitation rate. Therefore, we predict that the waves 
grow to nonlinear amplitudes, and the cosmic rays are scattered at nearly the maximum 
possible rate. The limits on B obtained under the assumption of scattering at the maximum 
rate are probably close to correct. 

Cosmic rays will interact dynamically with the thermal gas as long as they are strongly 
scattered [see eqns. (9) and (10)]. The rate of energy transfer is usually assumed to be 
proportional to the Alfven speed under standard conditions, but is much larger in the present 
case because the turbulence which scatters the cosmic rays is highly super-Alfvenic [see eqn. 
(44)]. Cosmic rays could potentially play an important role in driving outflows from young 
galaxies. 

The fact remains that magnetic fields and cosmic rays are in approximate equipartition 
in our own Galaxy at the present time. We suggest that this is so because they share a 
common energy source: supernovae. Given enough time, magnetic fields, cosmic rays, and 
for that matter turbulent bulk motions will probably achieve a state near equipartition, as 
currently observed in the Galaxy. However, as we have shown in this paper, this is not 
required by anything fundamental to the physics of cosmic ray acceleration or confinement. 
Why the magnetic field and cosmic ray energy density have saturated at values close to 
the turbulent energy density, instead of increasing beyond it, is not completely clear. The 
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reason may be that if the magnetic field or cosmic ray energy density were appreciably larger, 
the vertical stratification of the interstellar medium would be unstable to buoyancy driven 
instabilities (Parker 1966), which would probably lead to separation of the thermal and 
nonthermal components of the interstellar medium, and possibly facilitiate magnetic field 
and cosmic ray escape. We have briefly considered the role of these buoyancy instabilities 
at early times, and find it to be ambiguous. If the ratio of cosmic ray pressure to thermal 
gas pressure is denoted by a, the effective polytropic exponents of the thermal gas and 
cosmic rays are ^ g and 7 cr , and magnetic pressure is negligible, then the system is unstable if 
7 9 — 1 + a(7 cr — 1) < (Zweibel & Kulsrud 1975). Thus, if the cosmic rays are well coupled 
by scattering and 7 cr = 4/3, they can stabilize the system, but if they diffuse rapidly, they 
are destabilizing. While a detailed investigation of this problem is beyond the scope of this 
paper, it is safe to say that if the instability occurs, its nonlinear development differs from 
that of the instability in a strong magnetic field. 

If cosmic rays can be accelerated and confined in the presence of magnetic fields which 
are several orders of magnitude weaker than the fields in contemporary galaxies then evi- 
dence for cosmic rays in young galaxies is not evidence for magnetic fields at the equipartition 
level. The fieldstrengths required are several orders of magnitude larger than the seed fields 
produced by mechanisms operating on large scales, but are representative of the fields ex- 
pected from a superposition of many small sources, such as plerion supernova remnants. The 
dominance of small scale fields predicted by some dynamo theories is inconsistent with the 
presence of cosmic rays, unless the field is strong enough that the cosmic ray gyroradius is 
less than the size of the system. Thus, the lower bounds on B implied by the presence of 
cosmic rays impose meaningful constraints on theories of the origin and evolution of galactic 
magnetic fields (§2). 

Finally, the drastic changes in the properties of hydromagnetic turbulence brought about 
by cosmic rays may be of some importance beyond the present application, although they 
fall outside its immediate domain. It has been pointed out elsewhere that cosmic rays can 
be a significant source of hydromagnetic fluctuations at the scale of their gyroradius, and 
that this may enhance the rate at which they are accelerated in shocks [Eichler (1985), 
Lucek & Bell (2000), Bell & Lucek (2001)]. Because the gyroradius scale falls between the 
global galactic scale and the resistive scale, fluctuations driven by cosmic rays may affect the 
operation of the galactic dynamo. The novel dispersion relation of these fluctuations will 
affect their nonlinear behavior, and will modify the properties of a turbulent cascade. These 
effects could operate in any environment in which the magnetic field is weak and energetic 
particles are present. 

This paper was begun at the conference "Making Light of Gravity" held in honor of Sir 
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Martin Rees. I am grateful to the organizers for the opportunity to speak there. Fabian 
Heitsch, J. Michael Shull, and especially an anonymous referee provided useful comments on 
the manuscript. Material support was provided by NSF Grant 0098701 to the University of 
Colorado. 



Appendix A. Dispersion Relation with Cosmic Rays 



The dispersion relation for plane electromagnetic waves in a plasma follows from the 
equation for the first order electric field E 1 



kx(kx E{) + —KE X = 0, 



(Al) 



where the response of the plasma is contained in the dielectric tensor K. Equation (Al) and 
general expressions for the components of K are given in texts such as Krall & Trivelpiece 
(1973). 

We consider circularly polarized waves propagating parallel to the ambient magnetic 
field zB . In this case, K yy = K xx and K yx = —K xy , and eqn. (Al) reduces to 



2JU21 



K xx i iK xy 



c 2 k 



0J Z 



(E lxT iE ly )=0. 



(A2) 



The ± and =F signs denote the sign of circular polarization. The dispersion relation follows 
by setting the term in square brackets equal to zero. 



The dielectric terms are 



K xx i iK xy 



p 2 v (1 — fl 2 ) 
o — kvji ± u a 



dfoa fkv_ 
dp \uj 



0q 



p dfx 



dpdji, (A3) 



where the summation index a represents particle species and the fo a are the zero order 
phase space distribution functions averaged over gyrophase. The u ca are the relativistic 
gyrofrequencies. 

We evaluate eqn. (A3) for a system with four populations of particles: cold (zero 
temperature) protons with number density rii drifting with speed — zvp, cold electrons with 
the same number density and drift speed as the ions, proton cosmic rays which are isotropic 
in velocity space and have number density n cr , and an isotropic population of cold electrons, 
also with number density n cr . This choice guarantees that to zero order the system is 
electrically neutral and there is no net current flow, consistent with the assumption of a 
uniform magnetic field. We have chosen to work in the frame in which the cosmic rays are 
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isotropic, and transform the resulting dispersion relation to the rest frame of the plasma, 
because that is the simplest way to do the calculation. 

Foote & Kulsrud (1979) considered the effect of finite temperature (but did not include 
cosmic rays). As we argue in §5.2, Landau damping is the most important thermal effect. 
It requires wave propagation at an angle to B , and is considered in Appendix B. 



The distribution functions of the cold, drifting species are 



n 



foa = "45 (P - PDa) S(fi + 1), 



(A4) 



where PDa = n^a^D, and a represents protons or electrons. The distribution function of the 
cold, stationary electrons is 

/oa = ^(p), (A5) 



and we write the distribution function of the cosmic rays as 



her = 2 n cr0(p), 



(A6) 



where <f>(p) is normalized to unity. Substituting eqns. (A4), (A5), and (A6) into eqn. (A3) 
and carrying out the integrations gives the dispersion relation 



c 2 k 2 



uj s 



Annie 2 , , . 
1 =F — {uj + kv D ) 



±- 



4nn cr e 2 



UJ 

1 V 



1 



m e (uj ce ± uj ± kvo 
in pi 



+ 



1 



1 



rrii (u cp ± uj ± kvo) 



± 



m e {oj ce ± u) 4 rripUJ, 



( 2 2\ # . 



rripLVcp 4 



Pip 2 (1 - n 2 



-l VTPi/p dp 
where V denotes the principal part of the integral, 



dpdfi 



Pi = 



k 



(A7) 



(A8) 



is the minimum momentum at which the resonance condition eqn. (1) can be satisfied, and 
we have dropped uj compared with the other terms in the denominator of the last integrals. 
The right hand side of eqn. (A7) is K xx ± iK xy . 

The "1" on the right hand side of eqn. (A7) represents the displacement current. The 
first square bracket represents the thermal plasma. Since the gyrofrequencies uj cp , u ce much 
exceed uj and kvo, we approximate the electron term as l/m e uj ce and the proton term as 
1/mpUJcp ^(uj + kv £>)/m p uj 2 p . The dominant terms cancel: uj ce m e = —uj cp m p . The cancellation 
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occurs because in the low frequency waves considered here, the electron and ion motion is 
nearly the E x B drift, which gives no net current. 

The second square bracket represents the cosmic rays and the cold electrons that cancel 
their current. The first cosmic ray integral represents the resonant cosmic rays. It is this 
term which leads to instability if the relative drift of the ions and thermal plasma exceeds 
t>Ai- The second term represents the nonresonant cosmic rays. If we revert to the cold plasma 
approximation by taking the limit p\/p^> 1, this term reduces to l/m p uj cp , which cancels the 
electron term. The condition pi/p ^> 1 is equivalent to the condition that the wavelength 
of the wave is much larger than the particle gyroradius. If this condition is not fulfilled, the 
particles do not E x B drift, and their current does not cancel the cold electron current. 
This can lead to a relatively large effect on the dispersion relation even when n cr /rii <C 1. 

We simplify the cosmic ray term somewhat by integrating the resonant term by parts 
and performing the integration over fi in the nonresonant term. The cosmic ray integrals in 
the second square bracket of eqn. (A7) then reduce to 

1 TP/Pi 



m pi 

2 UlpUJ C p jp 



p4>dp ± 



Pi V 

m v u C p 4 



(p 2 - pf) In 



1 ±p/pi 



=F Zppi 



dp 



dp. 



(A9) 



We then incorporate the electrons directly into the nonresonant cosmic ray term by replacing 

l/m e (uj ce ± a;) by -l/m p u cp . 

Finally, we drop the displacement current and multiply eqn. (A7) by uo 2 v\ { /c 2 . The 
resulting dispersion relation is 

Tier 

rii 

1 =F p/pi 



k 2 v 2 



Ai — + kv D ) + UJUJ cp — X 



m r~ i r pip 

— / p\p<pdp — 



pi 



(p 2 - p\) In 



1 ±p/pi 



T 2pp l ± - 1 - 
3pi. 



^10) 



If vd = and p/pi is small but finite for the bulk of the particles, eqn. (A10) reverts to the 
dispersion relation for hydromagnetic waves in a hot plasma (Foote & Kulsrud 1979), with 
a factor of n cr /ni multiplying the thermal correction term. 

The quantity in {} is the function Q r introduced in eqn. (25). It is clear from eqn. (A10) 
that (i = —(*. The behavior of ( is illustrated in Figure 1 for the normalized distribution 
function 

(au) 

which is close to the Galactic cosmic ray spectrum for p/po ^> 1 (the observed power law 
index is approximately 4.6 rather than 4). Equation (23) follows from eqn. (A10) by shifting 
back to the rest frame of the thermal plasma. 
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Fig. 1. — Absolute value of Re(() (solid curve) and Im(Q (dashed curve) vs pi for the 
distribution function given in eqn. (All). 

Appendix B. Landau Damping Rate 



We calculate the Landau damping rate using the dispersion relation derived by Foote 
& Kulsrud (1979), but replacing the thermal terms entering into the real part of oo with 
cosmic ray terms. The dispersion relation for small propagation angle 9 can be written as 



c 2 k 2 



K xx + iK. 



XI) 



c 2 k 2 



K xx iK X y 



+ i 7r 1 / 2 — I K 



kvi 



c 2 k 2 \ c 2 



tan^ 9 = 0, 



(Bl) 

where = (2ksT ' /m p ) 1 / 2 , and the components of the dielectric tensor are given in Appendix 
A. For 9 = 0, eqn. (Bl) reduces to the dispersion relation eqn. (A2). 

We write uj = ojq + eoj\, where u)q is a solution of eqn. (Bl) for 9 = and assume 
e ~ 0(9 2 ) 1. Expanding eqn. (Bl) to first order in e and using the dispersion relation at 
9 = gives 



.7T 1 / 2 kVj C 2 



eui = —i- 



,2 
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V(u G y 1 tim 2 9 



where 



d 



D(uj ) = -q^ (Kxx ± iK xy ) Uo- 



(B2) 
(B3) 
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If the waves are Alfven waves, u 2 = k 2 v\ {) D = 2c 2 k 2 /uj 3 , and 

.7T 1 / 2 kvi , „ j . 

ecui=«— -tan 2 #, B4 

2 c<j 

which to order unity is the result of Foote & Kulsrud (1979). If the parameter X defined in 
eqn. (27) is large, to ~ (e/ r cj cp /a>£>) 1//2 [see eqn. (31)], and D ~ 2c 2 /u;'y 2 ii . Substituting this 
result into eqn. (B2) shows that eui is still given by eqn. (B4). 

Foote & Kulsrud (1979) found that the waves are critically damped - the real and 
imaginary parts of to are comparable - when 9 > (v Ai/^i) 1 ^ 2 - When the wave frequency is 
determined by cosmic rays, the angle for critical damping is increased to ~ (cv D^ir/v 2 ) 1 ^ ~ 
10(vi)/c) 1 / 4 (njT)~ 1 / 4 . The damping is weakened because the thermal plasma energy is a 
relatively small part of the wave energy, as shown in Appendix C. 



Appendix C. Damping Rate from Ion-Neutral Friction 

We calculate the damping rate Ti n using an energy method, rather than directly solving 
the dispersion relation Let the energy density in waves be U w and the rate per unit volume 
at which energy is dissipated be U w . These quantities are related to T in by 

U w = —2Ti n U w . (CI) 

The dissipation rate can be written in terms of the first order thermal ion and neutral 
velocities v Vl and v ln as 

U w = -Pifin | Vu - Vi n | 2 • (C2) 

The wave energy density includes contributions from the electric and magnetic fields, thermal 
plasma, and cosmic rays (the contribution from the neutrals is much smaller than the plasma 
contribution because the wave frequency is so high that the neutrals remain nearly at rest), 
and can be written as (Bekefi 1966) 

U w = ^-B\ ■ B 1 + ^-E\ ■ A (uK) ■ Ei, (C3) 

87T 87T OUJ 

where K is the dielectric tensor (see Appendix A). The first term in eqn. (C3) represents 
magnetic energy density, and the second accounts for the energy density in the electric field, 
thermal plasma, and cosmic rays. 

We assume the waves are circularly polarized and propagate parallel to B , which 
we take to define the z direction, and that first order quantities vary with z and t as 
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expi(ut — kz). The neutral velocity is determined solely by friction with the ions, and 
can be written in terms of Vu as 

iuv u = v ni (v ln - vu) , (C4) 

where the neutral- ion collision frequency u ni = p i (av)/(m i + m n ). The ion velocity is deter- 
mined by both frictional and magnetic forces (it suffices to use a one fluid treatment of the 
plasma, Vu ~ v± e , because the thermal electron and ion velocities are essentially the B x B 
drift. The ion equation of motion is 

iuvu = -- kB B 1 + v in (v ln - vu) . (C5) 

Solving eqn. (C4) for V\ n in terms of Vu and expressing Vu in terms of B±, we rewrite eqn. 
(C2) as 

U w = -u m V -^^B\-Bu (C6) 

where we have assumed uj /u in ^> 1 . 

The second term on the right hand side of eqn. (C3) can be written as 

±-E lx E* x -^-u (K xx ± iK xy ) . (C7) 

From eqns. (A2) and (31), 

v 2 Ai UJ 2 

It is straightforward to evaluate eqn. (C7) using eqn. (C8). We replace E\ X E\ X by B\ • 
Blv 2 /2c 2 using Faraday's law. Equation (C3) becomes 

«„ = ^ Bl . B: (l + |). (C9) 

The first term in eqn. (C9) represents the magnetic energy density, while the second is 
primarily cosmic ray energy, with a small contribution from the first order electric field. 

Equation (43) follows from eqn. (CI) if we substitute for U w and U w from eqns. (C6) 
and (C9). 
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